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SUMMARY 


Personnel  contributing  to  this  contract  are:  (1)  Dr.  Eugene  Herrin,  Principal 
Investigator,  (2)  Paul  Golden,  Project  Director,  (3)  Karl  Thomason,  Chief 
Engineer,  (4)  Nancy  Cunningham,  Director  -  Computer  Laboratory,  (5)  David 
Anderson,  Systems  Analyst,  (6)  Dyann  Anderson,  Administration,  (7)  Dick 
Arnett,  Consultant,  (8)  Herbert  Robertson,  Consultant,  (9)  Jack  Swanson, 
Consultant  (10)  Victoria  L.  Hammill,  Consultant  and  (11)  Billie  Myers, 
Consultant.  Ph.  D.  students  include:  (1)  Chris  Hayward,  (2)  Relu  Burlacu,  and 
(3)  Zenglin  Cui. 

Objectives 

Objectives  of  the  contract  are  twofold:  (1)  to  conduct  research  in  seismic  mini¬ 
array  technology  and  use  of  data  from  single  stations  and  sparse  networks, 
and  (2)  to  design,  evaluate,  and  construct  two  mini-arrays,  TEXESS  (Texas 
Experimental  Seismic  System)  in  Southwest  Texas  and  LUXESS  (Luxor 
Experimental  Seismic  System),  which  is  northeast  of  Luxor,  Egypt.  These  two 
tasks  are  dubbed  CLIN  1  and  CLIN  2. 

The  original  CLIN  1  objectives  were  to:  (1)  conduct  research  in  the  use  of 
single  station  and  sparse  network  data  in  detecting  and  identifying  small 
seismic  events,  (2)  conduct  research  to  develop  optimum  configurations  and 
processing  techniques  for  a  nine-element  mini-array,  and  (3)  to  continue 
development  of  an  unmanned  intelligent  seismic  station.  These  objectives 
have  been  revised  by  the  Project  Office  in  April  1994  as  described  on  page  6 
under  Implications  for  Further  Research. 

CLIN  2  objectives  are  to:  (1)  acquire  hardware  and  software,  (2)  install  TEXESS, 
(3)  perform  site  surveys  and  choose  location  for  LUXESS,  (4)  test  TEXESS  and 
perform  verification  tests  prior  to  de-installation,  (5)  de-install  TEXESS,  (6) 
complete  civil  work  in  Egypt,  (7)  install  and  test  LUXESS,  (8)  de-install  data 
acquisition,  analysis  and  archiving  equipment  and  ship  to  Helwan,  Egypt, 
data  center,  and  (9)  install  and  test  data  acquisition,  analysis  and  archiving 
equipment  at  Helwan  data  center. 
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Technical  Problem 


The  German  Experimental  Seismic  System  (GERESS)  was  dedicated  in  1992 
and  represents  an  upgrade  for  regional  arrays.  Although  GERESS  was 
technologically  advanced  over  NORESS  and  ARCESS,  which  were  earlier 
regional  arrays,  because  of  greater  sensitivity  and  wider  dynamic  range,  there 
was  a  considerable  effort  that  resulted  in  increased  costs  for  pier  and  vault 
construction  and  trenching  for  power  cabling.  Now  in  TEXESS,  there  are 
innovations  in  emplacement  techniques,  which  include  the  installation  of 
sensors  in  shallow  boreholes  instead  of  vaults  and  the  use  of  solar  power  at 
each  site  to  eliminate  cabling  from  a  central-power  source,  that  have  reduced 
array-installation  costs  by  an  order  of  magnitude.  TEXESS  is,  therefore,  a 
proposed  design  for  a  GSE-Alpha  station  because  of  these  cost-cutting 
innovations.  In  addition  to  design,  construction,  installation,  and  operation, 
of  TEXESS,  research  will  be  undertaken  to  develop  new  means  of  taking  data 
and  handling  the  data. 

General  Methodology 

In  GSE/US/84,  February  1993,  entitled  "Technical  Concepts  for  an 
International  Data  Exchange  System,"  the  GSE  established  the  design  goals  of 
a  future  system.  Goals  are: 

1.  Provide  prompt  access  to  all  essential  data, 

2.  Provide  convenient  access  to  all  available  data, 

3.  Provide  direct  access  to  all  data  at  authorized  national  and  global  facilities, 
and 

4.  Accomplish  goals  with  realistic  manpower  and  budget  resources. 

The  new  concept  of  a  global  system  for  data  exchange  calls  for  an  Alpha 
Network  of  40-60  stations,  primarily  arrays;  plus  much  greater  than  60 
Regional  or  Beta  Stations;  plus  Local  and  National  Networks  or  Gamma 
Stations. 
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SMU  began  research  on  mini-array  technology  in  1991  on  a  previous  contract. 
The  proposed  design  was  along  the  lines  of  an  Alpha  Station  consisting  of  an 
array  containing  nine  sites.  Advancements  over  the  GERESS  design  included 
the  following; 

1.  The  placement  of  seismometers  and  electronics  in  boreholes  to  greatly 
reduce  construction  costs  for  piers  and  vaults 

2.  The  use  of  solar  power  at  each  site  rather  than  a  central-power  source 

3.  The  use  of  GPS  receivers  for  time  data  at  each  seismometer  site  to  replace 
central  timing  from  the  Hub 

4.  The  employment  of  radio  links  from  seismometer  sites  to  the  Hub  to 
replace  cable  links  and  associated  construction  costs 

5.  The  use  of  modular  equipment  to  facilitate  the  installation  and 
maintenance  of  the  array. 

Four  shallow  boreholes  about  7  meters  deep  and  11 -5 /8-in.  in  diameter  were 
drilled  and  cased  with  standard  8-in.  pipe.  Special  equipment  and  techniques 
were  developed  to  lower  and  level  seismometers  in  the  boreholes.  A 
prototype  solar  power  array  and  directional  antenna  were  also  developed  for 
installation  at  LTX. 

Technical  Results 

The  limited  program  described  above  was  successful  and  SMU  was  granted  on 
26  April  1993  Contract  No.  F19628-93-C-0057  to  design,  evaluate,  and  construct 
two  nine  element  experimental  mini-arrays:  TEXESS  and  LUXESS. 

Important  Findings  and  Conclusions 

The  SMU  mini-array  research  program  that  was  begun  in  1991  under  the 
previous  contract  proved  the  feasibility  of  the  proposed  design  and 
methodology  described  above. 
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Significant  Hardware  Development 


Preliminary  research  has  led  to  the  following  hardware  developments: 

1.  The  development  of  seismometer  emplacement  techniques  in  boreholes, 
including  remote  seismometer  locking  eliminated  the  need  for  vaults 

2.  Advancements  in  computer  applications  and  radio  modems  allow  all 
necessary  electronic  components  to  fit  inside  a  8-in.  casing  to  provide  physical 
protection  and  a  more  stable  environment  for  the  electronics 

3.  The  use  of  Global  Positioning  Satellite  (GPS)  receivers  to  obtain  timing 
accurate  to  within  10  ps  of  world  time  assuring  time  synchronization  of  the 
array 

4.  The  use  of  modern  digital  radio  modems  allows  the  system  to  perform  as  a 
local  area  network  referred  to  as  a  RAN  (Radio  Area  Network);  radio  polling 
software  provides  wide  bandwidth  intra-array  communications  while 
requiring  two  base-station  radios;  the  need  for  expensive  buried  fiber-optic 
cable  is  eliminated 

5.  A  NEMA  enclosure  is  mounted  on  top  of  the  borehole  and  is  used  to  house 
the  batteries  and  as  a  mount  for  the  solar-power  array;  the  GPS  receiver  and 
radio  antenna  are  mounted  above  it. 

Special  Comments 

The  task  of  adapting  the  solar-panel  arrays  at  Lajitas  to  the  LUXOR 
environment  is  simplified  somewhat  in  that  both  TEXESS  and  LUXESS  are  at 
approximately  the  same  latitude,  30  deg  North;  both  are  in  arid  climatic 
zones;  and  both  have  about  3,500  annual  hours  of  sunshine.  As  a  result,  there 
would  be  no  need  to  modify  the  prototypic  TEXESS  design  because  of  differing 
environmental  conditions  at  LUXESS. 

Implications  for  Further  Research 

CLIN  1  objectives  were  revised  by  the  Project  Office  in  April  1994  to:  (1) 
conduct  research  to  develop  optimum  configurations  and  processing 
techniques  for  a  nine-  and  sixteen-element  short-period  arrays,  (2)  conduct 
research  in  discrimination  of  nuclear  events  using  autoregressive  (AR) 
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modeling  techniques  on  Lg  data,  and  (3)  conduct  research  in  measuring  20- 
second  Rayleigh  waves  at  regional  distances  using  high-resolution,  wide- 
dynamic-range,  short-period,  seismic-array  data  and  broadband  KS  36000  data. 

CLIN  1  -  RESEARCH 

Array  Research 

Conduct  research  to  develop  optimum  configurations  and  processing 
techniques  for  nine-and-sixteen  element  short-period  arrays, 

In  Scientific  Report  No.  1,  PL-TR-94-2106,  we  discussed  the  problems  of  the 
large  scatter  of  the  order  of  ±  15  deg  of  azimuth  estimates  at  GERESS  after  f-k 
processing.  In  order  to  address  this  problem,  SMU  research  has  concentrated 
on  developing  a  time-domain  processing  techniques  to  reduce  this  statistic 
using  the  nine-element  TEXESS  array  (see  Figure  3,  TEXESS  Map,  in  cited 
report).  The  array-processing  technique  is  similar  to  that  described  by  Bernard 
Massinon  in  his  paper  entitled  "The  French  seismic  network  -  current  status 
and  future  prospects,"  which  he  presented  at  the  GERESS  Dedication  and 
Symposium  on  24  June  1992.  The  processing  algorithm  developed  by  SMU 
using  GERESS  D-ring  data,  which  approximates  the  proposed  9-element 
TEXESS  array,  was  presented  in  SMU-R-92-396,  p.  14-17. 

Appendix  1  is  a  report  on  array  research. 

Discrimination  Research 

Conduct  research  in  discrimination  of  nuclear  events  using  autoregressive 
(AR)  modeling  techniques  on  Lg  data 

In  the  framework  of  a  Comprehensive  Test  Ban  Treaty  (CTBT), 
discrimination  between  low-yield  or  decoupled  nuclear  explosions,  economic 
explosions  and  small  shallow  earthquakes  using  the  characteristics  of  the 
seismic  waves  becomes  very  important.  Some  of  the  economic  explosions  are 
multiple-source  events  with  a  time  and  space  pattern  dependent  upon  the 
type  of  application.  The  superposition  of  the  seismic  motion  in  the  time 
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domain  leads  to  regular  amplification  and  suppression  of  spectral  power  in 
the  frequency  domain.  As,  in  general,  single  events  (single  explosions  or 
earthquakes)  do  not  exhibit  spectral  modulations,  their  presence  can  be  used 
in  the  discrimination  between  single  and  multiple  events.  The  aim  of  the 
present  study  is  to  develop  a  fast  and  robust  method  of  discriminating 
between  earthquakes  and  economic  explosions  based  on  differences  observed 
in  the  spectral  content  of  the  regional  waveforms.  The  method  is  based  on  the 
parametric  estimation  of  the  power-spectral  density  (PSD)  using  the 
autoregressive  (AR)  Burg  algorithm  of  order  3,  which  provides  a  fast  method 
to  emphasize  the  spectral  differences. 

Appendix  2  is  a  report  on  this  work. 

Conduct  Ms'.nih  research  by  measuring  20-second  Rayleigh  waves  at  regional 
distances  using  high-resolution,  wide-dynamic-range,  short-period,  seismic- 
array  data  and  broadband  KS  54000  data. 

The  Ms:mb  discriminant  has  been  investigated  by  a  number  of  researchers  for 
both  regional  and  teleseismic  events  and  explosions.  Bases  for  the 
discriminant  are  (1)  that  explosions  emit  more  energy  in  the  form  of  high- 
frequency  body  waves  and  (2)  that  earthquakes  emit  more  energy  in  surface 
waves  having  low  frequency  radiation;  therefore,  an  Ms:mb  plot  displays  a 
significant  separation  of  the  two  populations.  The  problem  with  the  method 
is  that  of  identifying  small  explosions;  that  is,  the  problem  boils  down  to 
seismograph  sensitivity.  With  the  installation  of  new  high-dynamic-range 
seismographs  at  TEXESS,  planned  research  includes  the  determination  of  Ms 
from  small  earthquakes  at  regional  distances  using  the  TEXESS  array  data 
recorded  by  short-period  GS-13  seismometers  and  a  posthole,  long-period  KS 
54000  seismometer. 

Victoria  L.  Hammill  has  begun  development  of  discrimination  software  to 
provide  an  automated  system  of  analysis  for  event  discrimination  through 
the  development  of  public  domain  seismic  applications  software,  with 
particular  emphasis  on  conversion  of  cumbersome  CALPLOT  graphics  to 
more  efficient  Matlab  graphics. 
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Appendix  3  is  a  report  on  Ms:mb  research,  and  Appendix  4  is  a  Poster  Session 
on  discrimination  research  at  SMU. 

CLIN  2  -  DESIGN,  EVALUATION,  AND  CONSTRUCTION  OF  TEXESS  AND 

LUXESS 


Experimental-Array  Program 

Information  on  the  experimental-array  program  at  SMU  on  the  previous 
contract  was  presented  in  SMU-R-92-396  and  in  Scientific  Report  No.  1,  PL- 
TR-94-2106  (ADA284580). 

TEXESS  AND  LUXESS 

Acquisition  of  Hardware  and  Software 

The  First  and  Second  Quarterly  R  &  D  Status  Reports  cover  the  acquisition  of 
hardware  and  software.  TEXESS  and  LUXESS  equipment  is  discussed  in 
Scientific  Report  No.  1,  PL-TR-94-2106.  Instructions  for  the  installation  of  the 
Posthole  54000  seismometer  are  presented  in  Appendix  5  of  this  report. 

Array  Hardware 

Hardware  is  discussed  in  Scientific  Report  No.  1,  PL-TR-94-2106. 

Computer  Hardware 

Computer  equipment  is  discussed  in  Scientific  Report  No.  1,  PL-TR-94-2106. 
Software 


Acquisition  of  software  was  addressed  in  Scientific  Report  No.  1,  PL-TR-94- 
2106. 
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Install  TEXESS 


Layout 

TEXESS  layout  is  discussed  in  Scientific  Report  No.  1,  PL-TR-94-2106. 
Installation 


Installation  is  discussed  in  Scientific  Report  No.  1,  PL-TR-94-2106. 

Perform  Site  Surveys  and  Choose  Locations  for  LUXESS 

SMU  has  received  communications  from  the  National  Research  Institute  of 
Astronomy  and  Geophysics  (NRIAG),  Helwan,  Egypt,  requesting  a  draft 
Memorandum  of  Understanding  (MOU).  Arrangements  are  presently  being 
pursued  by  respective  State  Departments. 

Test  TEXESS  Prior  To  De-installation 

System  fidelity  tests  were  conducted  by  Chris  Hayward  and  Dick  Kromer 
during  the  week  of  14  November  1993.  During  this  period,  engineering 
refinements  were  made  by  Karl  Thomason. 

De-install  TEXESS 

During  the  early  part  of  the  week  of  21  November,  SMU  principals  made  an 
inspection  trip  to  TEXESS  in  order  to  make  plans  for  de-installation,  packing, 
and  international  shipment. 
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APPENDIX  1.  -  ARRAY-PROCESSING  RESEARCH 


A  flow  diagram  for  the  array-processing  algorithm  is  presented  in  Figure  1. 
After  computing  the  cross-correlations  of  all  channels  relative  to  a  reference 
channel  and  determining  the  indices  of  maximum  correlation  for  each 
channel,  relative  time  delays,  At's,  are  computed.  Using  trial  azimuths  and 
incrementing  by  1  degree,  time  residuals  are  computed  and  the  rms  residual 
is  calculated  for  each  azimuth.  The  process  is  repeated  using  each  of  the  nine 
stations  as  reference  channels.  The  trial  azimuth  yielding  the  minimum  rms 
residual  is  selected  as  the  true  azimuth.  Figure  2  shows  the  standard 
deviation,  SD,  of  the  true  azimuth  is  1.4  deg  when  using  all  nine  channels  as 
reference. 

Having  selected  the  best  azimuth,  one  station  is  chosen  as  the  reference 
channel  and  the  time  differences  are  corrected  for  elevation,  using  a  trial 
replacement  velocity.  The  differences  are  then  regressed  against  the  distances 
from  the  reference  channel  and  the  slope  of  the  line  gives  the  velocity,  Vh, 
and  the  residuals  for  each  station  are  then  computed  and  regressed  against  the 
elevation  differences  determining  the  best  replacement  velocity,  Vr,  as  shown 
in  Figure  3.  This  procedure  is  then  repeated  using  each  station  as  the 
reference  channel  and  the  velocities,  Vh  and  Vr,  are  averaged  to  yield  the 
true  velocities  to  aid  in  the  determining  the  phase  of  the  arrival  as  shown  in 
Figure  4.  Figure  5  shows  standard  deviations  of  azimuths  of  7.1  obtained  by  /- 
k  processing  of  the  25-element  GERESS,  which  decreases  to  5  deg  using  time- 
domain  processing  of  9-elements  of  the  GERESS  array  to  simulate  the  TEXESS 
geometry. 

For  the  9-element  TEXESS,  cross-correlation  techniques  in  the  time  domain 
were  used  to  determine  relative  time  delays,  At's,  between  array  elements  for 
an  event.  After  obtaining  the  At's,  an  iterative  procedure  was  used  to 
determine  the  azimuth,  station  correction,  and  phase  velocity  for  the  event. 
Time-domain  processing  eliminates  the  dual  problems  of  space  aliasing  and 
side  lobes  of  the  main  beam  that  are  inherent  in  frequency-domain 
processing. 
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Initialize  velocities 
V  =8.50km/s  V  =  6.50  km /s 


Loop  until  velocities  stabilize 


Initialize  ref  channel  to  1 
Initialize  vel's  to  current  estimates 


Loop  over  channels  as  ref^ 


Residual  rms  (secs) 


Residual  rms  vs  trial  azimuth,  1992148_1814_16 


Azimuth  (deg) 


Figure  2.  —  Trial  Azimuth  vs  Residual  rms  Using  9-element  GERESS  D-ring. 

1  1 


Time  difference  corrected  for  elevation  (secs)  Time  difference  corrected  for  distance  (secs) 


Pn  and  replacement  velocity  estimates,  Pack  azimuth  =  194  peg 


Figure  3.  —  Regression  Lines  to  Obtain  Vh  and  Vr. 


Figure  4.  —  Regression  Lines  to  Obtain  True  Vh  and  Vr. 


Number  of  events  Number  of  events 


Deviation  of  25-element  GERESS  azimuth  from  PDE  azimuth.  1992  events 


Deviation  of  9-element  GERESS  azimuth  from  PDE  azimuth,  1992  events 


Figure  5.  -  Histogram  of  Time  Residuals  Showing  Decrease  in  Standard 
Deviation  of  Azimuths  Using  Time-domain  Method  Compared  with  f-k 

Method. 


The  following  illustrations  show  the  above  technique  applied  to  a  magnitude 
5.0  aftershock  recorded  by  TEXESS  on  19  January  of  the  mb=6.8  Northridge, 
California,  earthquake  of  17  January  1994.  Figure  6  shows  the  P-wave  arrivals 
bandpass  filtered  from  0.5  to  5.0  Hz.  Time  delays,  At's,  are  computed  from 
cross-correlations  in  a  2-sec-long  time  window.  Figure  7  shows  the  rms  error 
as  a  function  of  assumed  azimuth  for  reference  site.  Cl,  using  a  P-wave 
velocity  of  8.43  km/ sec.  The  minimum  rms  error  corresponds  to  an  azimuth 
of  294°.  Figure  8  shows  the  beam  sum  of  the  8  vertical  elements  of  the  array 
with  the  beam  steered  at  294°  using  a  velocity  of  8.429  km/ sec.  Figure  9  shows 
the  P-wave  velocity  determined  by  a  least-squares  fit  to  the  residuals  vs 
distance  differences  relative  to  the  reference  site.  Cl,  using  an  azimuth  of 
294°.  Figure  10  is  the  distribution  of  the  time  residuals  using  the  same 
azimuth  and  a  velocity  of  8.43  km /sec. 


C5 

C4 

C3 

C2 

Cl 

B3 

B2 

B1 


AO 


1994019  20_00  mb  5.0  Pn  0.5-10  Hz 


Figure  6.  -  TEXESS  P-wave  Arrival  for  mb=5  Northridge,  CA,  Aftershock  of 
19  January  1994  Eollowing  the  Northridge  mb=6.8  Quake  of  17  January  1994. 


RMS  error  (secs) 


1994019  20_00  mb  5.0  Pn  rms  error  vs  azimuth  reference  site  Cl 


Figure  7.  —  Azimuth  vs  rms  Error. 
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ure  8.  -  TEXESS  Beam  Sum. 


Fig 


Elevation-corrected  residual  (secs) 


1994019  20_00  mb  5.0  P  least-squares  fit  to  residuals  vs  distance  differences  reference  site  Cl 


Figure  9.  —  Distance  Differences  vs  Residuals. 


APPENDIX  2.  --  SPECTRAL  DISCRIMINATION  USING  A  PARAMETRIC 

APPROACH 


Introduction 

Spectral  content  of  regional  phases  (Pn,  Pg,  Lg)  has  been  shown  to  be 
important  in  the  design  of  algorithms  for  discrimination  between 
earthquakes  and  explosions  (Murphy  and  Bennett,  1982).  Techniques  using 
spectral  ratios  have  been  applied  to  data  in  the  Western  U.S.  and  in  Europe 
with  reasonable  success,  but  the  transportability  of  the  methods  has  not  been 
satisfactorily  established.  Wuster  (1993)  used  ARMA  (Auto-Regressive, 
Moving  Average)  spectral  estimates  in  the  development  of  discriminants. 
All  spectral  discriminants  based  on  regional  arrivals  have  at  this  time  at  least 
two  deficiencies: 

(1)  They  have  not  been  adequately  tested  against  data  with  a  variety  of 
sources  and  propagation  paths,  and 

(2)  There  is  not  an  adequate  explanation  for  why  they  work  based  on 
empirical  or  theoretical  source  models. 

Research  is  underway  to  eliminate  these  deficiencies  at  least  for  one  spectral 
discriminant  that  has  been  shown  to  work  very  well  in  central  Europe. 

Background 

The  discriminant  we  propose  to  study  is  based  on  autoregressive  modeling  of 
the  Lg  phase.  The  autoregressive  spectral  model  is 


x(n)  =  -^a{k)x{n-k)  +  u{n) 

where  p  is  the  order  of  the  process,  a(k)  are  the  coefficients 
and  u(n)  is  the  white  noise  driving  process. 

In  the  study  discussed  here  AR  order  3  was  used.  The  resulting  spectral 
model  consists  of  one  complex  paired  pole  and  one  real  pole,  the  latter  being 
located  at  zero  frequency  for  the  events  studied  so  far.  (A  real  pole  can  only  be 
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located  at  zero  frequency  or  at  the  folding  frequency.)  The  initial  data  is  from 
the  CSS  Ground  Truth  Data  Base  (Grant  et  al.,  1993)  recorded  at  GERESS  (16 
Hz  bandwidth).  The  data  consists  of  explosions  and  earthquakes  in  the 
Vogtland  area  which  are  tabulated  in  Table  1.  Figure  1  shows  the  third  order 
autoregressive  spectrum  (AR(3))  for  a  quarry  blast  and  Figure  2  shows  the  AR 
(3)  spectrum  for  an  earthquake.  The  frequency  and  reciprocal  pole  position  of 
the  complex  pole  in  the  AR  (3)  models  were  calculated  using  the  Lg  arrival 
for  the  events  in  Table  1.  A  value  of  1.0  indicates  that  the  reciprocal  pole  is  on 
the  unit  circle,  0.0  indicates  a  position  at  the  origin.  Neither  a  value  of  0.0  or 
1.0  can  occur  in  practice.  Figure  3  shows  reciprocal  pole  position  versus 
frequency  for  the  analyzed  events.  Solid  circles  are  explosions  and  stars  are 
earthquakes.  The  two  open  circles  represent  surface  and  buried  explosions  at 
the  Tyrnyaus  mine  in  the  Caucasus  Mountains  of  Southern  Russia  (Stump  et 
al.,  1994)  recorded  at  station  GUM.  This  figure  shows  a  clear  separation  of 
explosions  and  earthquakes  with  the  latter  having  broad  AR(3)  spectra  with 
“weak"  poles  above  6  Hz  while  the  explosions  all  show  much  “stronger" 
poles  at  frequencies  less  than  5  Hz.  After  reviewing  the  results  shown  in 
Figure  3,  David  Russell  of  AFTAC  made  the  following  comment  (electronic 
communication,  1994): 

I  must  tell  you  I  am  very  impressed  with  this  approach,  since: 

(1)  It  appears  to  work 

(2)  It  is  robust,  simple  and  straightforward  (music  to  our  ears 

here) 

(3)  Last  but  certainly  not  least;  it  does  not  require  extensive 
precalibration, 

such  as  taking  out  Q  effects,  etc. 


Unresolved  Technical  Issues 

The  Lg  discriminant  based  on  a  third  order  autoregressive  model  has  been 
successfully  applied  to  the  available  seismic  data  with  a  bandwidth  of  at  least 
15  Hz.,  where  the  nature  of  the  source  was  clearly  identified.  Technical  issues 
that  now  must  be  resolved  can  be  simply  described  as  follows: 

(1)  Will  the  method  work  for  sources  in  other  areas,  particularly  areas 
analogous  to  and  within  the  Middle  East? 
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Table  1.  Earthquakes  and  Explosions  From  Vogtland  Area 
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Figure  2.  --  Seismogram  ar\d  PSD  AR(3)  for  EVENT  6. 
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Figure  3.  --  Absolute  Pole  Position  as  a  Function  of  Frequency. 
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(2)  The  method  has  been  tested  only  for  events  with  ML  around  2. 
Will  the  discriminant  work  equally  well  for  larger  events? 

(3)  Why  does  the  method  work  at  all?  That  is,  why  does  the  vertical 
component  of  the  Lg  signals  recorded  on  systems  having  a  flat  velocity 
response  from  about  1  to  15  Hz  show  an  apparent  "corner  frequency"  of 
less  than  5  Hz  for  shallow  distributed  explosions  and  greater  than  6  Hz 
for  earthquakes? 
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APPENDIX  3.  -  DETERMINATION  OF  MS  FROM  REGIONAL  EVENTS 


Introduction 

By  1970,  research  in  the  seismic  discrimination  of  explosions  and  earthquakes 
had  clearly  established  that  the  comparison  of  surface  wave  magnitude  (Ms) 
and  body  wave  magnitude  (mb)  could  be  used  to  classify  shallow  events.  In 
fact,  this  method  provided  the  only  reliable  means  for  the  seismic 
identification  of  underground  explosions.  The  Ms:mb  method  was  shown  to 
be  an  effective  and  transportable  discriminant  for  shallow  events  with  mb 
greater  than  4.75  (Davies,  1968).  Further  research  was  carried  out  including 
the  construction  of  long-period  arrays  in  an  attempt  to  lower  this  threshold, 
but  until  recently,  essentially  all  data  was  obtained  from  stations  at  teleseismic 
distances. 

Now  the  construction  of  regional  arrays  for  monitoring  a  proposed  CTBT  has 
led  to  considerable  interest  in  reducing  the  detection  threshold  for 
fundamental  mode  Rayleigh  waves  using  regional  signals.  The  goal  of  our 
proposed  studies  is  to  identify  shallow  earthquakes  at  regional  distances  with 
magnitude  (mb)  as  low  as  3.0  using  the  Ms:mb  method. 

Previous  Work 

The  GERESS  regional  array  in  eastern  Bavaria  provided  the  first  high- 
resolution  array  data  with  linear  dynamic  range  greater  than  120  dB. 
Although  in  1992  there  was  no  broadband  sensor  at  GERESS  which  could 
resolve  long-period  noise  at  a  20  sec.  period,  it  was  possible  to  use  the  short- 
period  sensors  for  this  purpose.  The  GS-13  sensors  at  GERESS  have  their 
resolution  at  0.05  Hz  (20  sec.  period)  limited  by  self-noise  which  is  about  20  dB 
above  the  Low  Noise  Model  (Peterson,  1993,  and  Rodgers,  1992).  The  random 
nature  of  this  self-noise  and  the  small  aperture  of  the  array  relative  to  the 
wavelength  of  20  sec.  Rayleigh  waves  make  it  possible  to  improve  the  signal- 
to-noise  ratio  by  about  14  dB  simply  by  summing  the  output  of  the  25  short- 
period  sensors.  The  120  dB  resolution  then  allows  the  sum  to  be  low-pass 
filtered  (corner  at  0.1  Hz)  in  order  to  observe  long-period  signals.  This  method 
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has  been  used  to  resolve  the  ambient  seismic  background  noise  at  GERESS  at 
periods  of  10  to  20  seconds. 

The  prototype  AREA  Model  94  regional  array  at  Lajitas  (TEXESS)  has  nine 
vertical  GS-13  (SP)  elements  and  a  post-hole  KS36000  broadband  (BB)  system. 
The  SP  and  BB  data  are  recorded  with  linear  dynamic  range  greater  than  120 
dB,  with  the  SP's  using  SHI  AIM  24  ADCs  and  the  BB  using  an  RDAS  200 
ADC.  The  prototype  post-hole  KS36000  at  TEXESS  is  sandpacked  in  a  7  m 
deep  borehole  resulting  in  a  total  installed  cost  less  than  half  that  of  the 
standard  KS36000  in  a  100  m  deep  borehole.  An  AFT  AC  standard  KS36000 
system  is  collocated  with  the  post-hole  KS36000,  but  the  data  from  that 
instrument  are  recorded  with  a  linear  dynamic  range  of  about  70  to  80  dB. 

The  sum  of  the  GS-13  elements  at  TEXESS  results  in  a  self-noise  threshold 
about  10  dB  above  the  Low  Noise  Model  at  a  period  of  20  sec.  The  posthole  KS 
should  be  capable  of  resolving  the  Low  Noise  level  at  20  sec.  We  have  used 
both  the  SP  array  sum  and  the  posthole  BB  vertical  to  look  for  20  sec.  Rayleigh 
waves  in  regional  signals. 

Aftershocks  from  the  Northridge  earthquake,  a  distance  just  over  1500  km 
from  Lajitas,  were  used  to  prepare  a  database  for  which  location,  origin  time, 
local  magnitude  and,  in  some  cases,  mb  were  known  from  USGS  reports. 
One  aftershock  on  19  January  1994  (mb  5.0)  was  selected  as  a  reference  event. 
First  the  SP  array  was  beamed  at  the  event  using  a  cross-correlation  method  to 
determine  back  azimuth,  horizontal  phase  velocity  and  static  corrections  for 
the  Pn  arrival  (az.  291’,  vel.  8.46  km/sec.).  Figure  1  shows  that  the  residuals 
for  this  fit  for  all  possible  cross-correlations  between  elements  are  all  less  than 
one  sample  point  (0.025  sec.  for  40  sps  data).  The  set  of  delays  from  this  fit  was 
used  to  form  beams  for  24  smaller  aftershocks.  These  array  sums  were  then 
used  to  compute  mb  (Pn)  for  Lajitas,  which  are  shown  in  the  eighth  column 
of  Table  1. 

Both  the  sum  of  the  SP  channels  and  the  vertical  BB  (post-hole)  channel  were 
Butterworth,  low-pass  filtered  at  0.10  Hz  with  three  poles  forward  and  three 
reversed  in  order  to  eliminate  any  phase  shift.  The  filtered  wave  forms  for 
the  SP  sum  and  BB  Z  were  essentially  identical.  Figure  2  shows  the  results  of 
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Number 


1994019  20  00  mb  5.0  P  distribution  of  residuals 


Residual  (secs) 


Histogram  of  residuals  from  least-sQuares  estimate  of  azimuth  and 
phase  velocity  of  reference  event  1994019_20_00  (aftershock  of  the 
Northridge  earthquake)  observed  on  the  Lajitas,  TX  short  period 
array.  Data  are  combined  for  residuals  obtained  using  each  channel  as 
reference.  Estimated  azimuth  was  291  degrees  and  estimated  phase 
velocity  was  8.46km/s. 


Figure  1.  —  Histogram  of  Residuals  from  Northridge  Aftershock. 
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Multiple-filter  plot  of  the  reference  event  sig94019_20_00,  nib=5.0 
ML=5.3,  showing  the  fundamental  Rayleigh  wave  group  velocity 
dispersion  curve.  This  event  was  used  to  develop  the  phase-matched 
filter  for  the  Northridge,  California  to  Lajitas,  Texas  travel  path. 

Wave-form  to  the  far  right  is  the  input  waveform  plotted  on  a  linear  time 
scale.  The  inside  waveform  is  the  input  signal  plotted  on  a  non-linear  time 
scale  corresponding  to  group  velocity  arrival  on  the  multple-filter  plot. 


Figure  2.  -  Rayleigh  Wave  Group  Dispersion  Curve  for  Northridge 

Aftershock. 
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a  multiple  filter  analysis  of  the  reference  event  using  a  program  supplied  by 
R.  Herrmann.  The  indicated  Rayleigh  wave  dispersion  curve  was  used  to 
design  an  optimum  phase-matched  filter  (Herrin  and  Goforth,  1977)  which 
when  applied  to  the  reference  event  produced  the  results  shown  in  Figure  3. 
This  filter,  adjusted  for  slight  differences  in  epicentral  distance,  was  applied  to 
the  low-pass  filtered  SP  sum  for  24  events  and  the  low-pass  filtered  BB  Z  for  6 
events.  Spectral  Ms  at  0.05  Hz  was  then  calculated  from  the  Pseudo 
Autocorrelation  Function  (PAF)  produced  by  the  phase-matched  filters  using 
the  formula 


Ms=  log  A/T  +  log  A  (1) 

where  A  is  the  spectral  amplitude  and  A  is  the  epicentral  distance  in  degrees. 

Table  1  shows  the  mb  (Pn)  and  Ms  values  for  the  25  events  which  are  ordered 
in  decreasing  values  of  ML.  It  was  observed  that,  on  the  average,  mb  was 
about  0.3  magnitude  units  smaller  than  ML  so  that  the  range  of  magnitudes 
for  the  data  base  was  from  about  mb  2.7  to  5.0. 

Ms  values  from  the  SP  sum  and  the  BB  Z  were  essentially  the  same.  For  both 
data  sources  the  20  sec.  Rayleigh  wave  was  lost  in  the  noise  for  events  with 
Ms  less  than  about  3.0  (mb  about  3.5).  From  Table  1  it  is  clear  that  Ms  is 
decreasing  faster  than  mb  with  decreasing  event  size.  We  can  speculate  that 
the  value  of  Ms  (20  sec.)  for  a  mb  3.0  Northridge  aftershock  would  be  about 
2.5,  but  neither  the  SP  sum  or  the  BB  Z  were  able  to  resolve  20  sec.  Rayleigh 
waves  at  this  level. 

Unresolved  Technical  Issues 

The  1500  km  path  from  Northridge  to  Lajitas  crosses  the  Coast  Ranges,  the 
Basin  and  Range  Province,  and  the  southward  extension  of  the  Rio  Grande 
Rift.  This  path,  roughly  analogous  to  paths  across  tectonically  active  areas  in 
the  Middle  East,  causes  considerable  attenuation  of  Pn  and  Rayleigh  waves. 
Thus,  the  use  of  the  Northridge  -  Lajitas  data  set  to  test  the  Ms:mb  technique 
at  regional  distances  constitutes  a  kind  of  “worst-case"  test  of  the  method. 
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Relative  Amplitude 


Reference  Event  siG94019_20_00  mb=5.0  ML=5.3 


Results  of  phase-matched  filtering  of  the  reference  event: 

Top  trace  is  the  input  waveform  of  reference  event  sig94019_20_00. 
The  middle  trace  shows  the  estimated  fundamental  Rayleigh  wave 
component  of  the  complex  input  signal  (pseudo  autocorrelation  function 
PAF).  The  bottom  trace  is  a  plot  of  the  residual  after  the  estimated 
fundamental  Rayleigh  wave  (middle  trace)  is  removed  from  the  input 
signal. 


Figure  3.  -  Phase  Matched  Filtering  of  Northridge  Aftershock. 
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If  we  are  to  lower  the  Ms  measurement  threshold  for  this  data  set  to  mb  3.0, 
we  must  investigate  the  nature  of  the  noise  floor  which  now  limits  the 
resolution.  The  SP  sum  at  20  sec.  must  have  a  noise  floor  at  or  greater  than 
the  instrument  noise  (lOdB  above  the  Low  Noise  Model)  and  the  BB  Z  must 
have  a  noise  floor  set  by  the  ambient  20  sec  noise  at  Lajitas  (near  the  Low 
Noise  level)  or  by  non-linear  effects  caused  by  the  large  output  voltages  at 
higher  frequencies.  But  the  SP  sum  and  BB  Z  were  recorded  on  different 
systems  and  both  systems  exhibited  the  same  noise  floor  well  above  the 
predicted  limits;  that  is,  essentially  equivalent  to  the  spectral  amplitudes  for 
Ms  3.0  events. 

What  is  the  nature  of  the  20  sec.  "noise"  which  limits  the  Ms  measurements 
at  Lajitas  for  the  Northridge  events?  How  is  the  "noise"  related  to  ambient 
background  in  the  absence  of  regional  signals?  Can  methods  be  found  for 
increasing  the  signal-to-noise  ratio  by  about  10  dB  in  order  to  obtain  Ms 
measurements  for  mb  3.0  events  in  Southern  California? 

For  Northridge  events  of  mb  about  3.5  the  Ms  (20  sec.)  value  has  decreased  to 
about  3.0.  The  use  of  phase-matched  filtering  has  removed  the  dispersion 
effects  so  that  spectral  amplitude  measurements  at  or  near  an  Airy  phase  can 
be  made  reliably.  Quoting  from  Davies  (1968,  SIPRI,  p.  62)  "When  magnitude 
determination  at  20  seconds  proves  impossible  at  near  distances,  Thirlaway 
considers  12  sec.  period  waves  and  applies  an  appropriate  correction...." 
Would  a  choice  of  12  or  15  sec.  period  for  the  Ms  measurement  result  in  a 
significantly  lower  threshold? 
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MsrnU)  STUDIES 


By  1970,  research  in  the  seismic  discrimination  of  explosions  and 
earthquakes  had  clearly  established  that  the  comparison  of  surface  wave 
magnitude  (Ms)  and  body-wave  magnitude  (mb)  could  be  used  to  classify 
shallow  events.  In  fact,  this  method  provided  the  only  reliable  means  for 
the  seismic  identification  of  underground  explosions.  The  Ms:mb  method 
was  shown  to  be  an  effective  and  transportable  discriminant  for  shallow 
events  with  mb  greater  than  4.75.  Further  research  was  carried  out 
including  the  construction  of  long-period  arrays  in  an  attempt  to  lower 
this  threshold,  but  until  recently,  essentially  all  data  was  obtained  from 
stations  at  teleseismic  distances. 

Now  the  construction  of  regional  arrays  for  monitoring  a  proposed  CTBT 
has  led  to  considerable  interest  in  reducing  the  detection  threshold  for 
fundamental  mode  Rayleigh  waves  using  regional  signals.  Our  goal  is  to 
identify  shallow  earthquakes  at  regional  distances  with  magnitude  (mb)  as 
low  as  3.0  using  the  Ms:mb  method. 
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1 99401 9_20_00  Pn  ML  5.3  1534  km  0.75  -  5.0  Hz 


Figure  I.  -  Pn  arrival  1994019_20_00  and  1994029_08_00 

Beam  of  Pn  arrival  for  two  aftershocks  of  the  Northridge  earthquake. 
Beams  were  formed  using  time  delays  determined  for  the  reference 
event  1994019  20  00.  Expected  Pn  arrival  times,  based  on  epicentral 
distances  reported  by  the  USGS  and  Pn  group  velocity  of  7.69  km/s, 
are  indicated  on  the  plots.  The  USGS  reported  mjj  5.0  and  no  ML  for 
1994019  20  00  and  m^  4.5  and  ML  5.0  for  event  1994029  08  00. 
ML  for  1994019  20  00  was  estimated  as  5.3,  based  on  the  relation 
developed  in  this  study. 
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1 99401 9_20_00  mb  5.0  P  distribution  of  residuals 

12  I - ^ ^ ^ ^ i 


Residual  (secs) 


Figure  2.  --  Histogram  of  residuals  from  least-squares  estimate  of 
azimuth  and  phase  velocity  of  reference  event  1994019_20_00 
(aftershock  of  the  Northridge  earthquake)  observed  on  the  Lajitas,  TX 
short  period  array.  Data  are  combined  for  residuals  obtained  using 
each  channel  as  reference.  Estimated  azimuth  vt^as  291  degrees  and 
estimated  phase  velocity  was  8.46km/s. 
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Figure  3.  --USGS  Magnitudes  for  Northridge  aftershocks 

Comparison  of  body  wave  and  local  magnitudes  reported  by  the 
USGS(GS)  for  8  aftershocks  of  the  Northridge  earthquake.  mb(GS) 
averages  0.3  magnitude  units  smaller  than  ML(GS). 


Comoarison  of  Lajitas  broadband  to  shod  period  20  second  spectral  amplituoes 


Short  period  ATT  (nm/s) 


Figure  4.  --  Comparison  of  broadband  to  short  period  20  second 
spectral  amplitudes  observed  at  Lajitas,  TX  for  7  aftershocks  of  the 
Northridge  earthquake.  Agreement  of  amplitudes  observed  on  the  two 
systems  is  good.  Gains  for  the  short  period  system  were  132.8 
counts/(nm  /s)  at  a  period  of  1  second  and  0.35  counts/(nm/s)  at  a 
period  of  20  seconds.  For  the  broadband  system,  gains  were  2.0E4 
counts/(nm/s)  at  a  period  of  1  second  and  466.8  counts/(nm/s)  at  a 
period  of  20  seconds. 


GROUP  VELOCITY  (KM/SEC) 


sic  94  019_  2.0_00sum_dl0 

DT  =0.2  5  0  SEC  DIST=15?4,16  NSAMP  =  2048 
c  AO  =50.270  A1  =0.000  AZ  =106.6 


PERIOD  (SEC) 


Figure  5.  --  Multiple-filter  plot  of  the  reference  event  sig94019_20_00, 
mb=5.0  ML=5.3,  shov>^ing  the  fundamental  Rayleigh  wave  group 
velocity  dispersion  curve.  This  event  was  used  to  develop  the  phase- 
matched 

filter  for  the  Northridge,  California  to  Lajitas,  Texas  travel  path. 

Wave-form  to  the  far  right  is  the  input  waveform  plotted  on  a  linear  time 
scale.  The  inside  waveform  is  the  input  signal  plotted  on  a  non-linear  time 
scale  corresponding  to  group  velocity  arrival  on  the  multple-filter  plot. 
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'A/VV 


Relative  Amplitude 


Reference  Event  sig94019_20_00  mb=5.0  ML=5.3 


Figure  6.  --  Results  of  phase-matched  filtering  of  the  reference  event: 

Top  trace  is  the  input  waveform  of  reference  event  sig94019_20_00. 
The  middle  trace  shows  the  estimated  fundamental  Rayleigh  wave 
component  of  the  complex  input  signal  (pseudo  autocorrelation  function 
PAF),  The  bottom  trace  is  a  plot  of  the  residual  after  the  estimated 
fundamental  Rayleigh  wave  (middle  trace)  is  removed  from  the  input 
signal. 
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LA  aftershock  mac  5.0  1 99401 9_20  summation  0.0-0. 1  Hz  winaowed  PAF 


SUMMARY 

Ms:mb  FOR  NORTHRIDGE  AFTERSHOCKS 


^  sp  array  and  bb  surface  wave  magnitudes  are  essentially  the  same. 

^  uses  mjj  is  on  the  average  about  0.3  units  less  than  ML  for  the 

Northridge  aftershocks. 

^  Ms  could  be  reliably  determined  for  events  with  mjj  about  3.5  and 
greater  at  a  distance  of  about  1500  km  from  Lajitas. 

^  Lajitas  m^,  was  determined  using  an  array  beam  with  an  azimuth  of 
272°  and  a  horizontal  phase  velocity  of  8.4  km/s. 

^  Pn  arrivals  for  Lajitas  with  m^j  less  than  about  3.3  were  heavily 
contaminated  by  background  noise. 


AR(3)  -  Lg  METHOD 


Spectral  content  of  regional  phases  (Pn,  Pg,  Lg)  has  been  shown  by 
Murphy  and  Bennett  to  be  important  in  the  design  of  algorithms  for 
discrimination  between  earthquakes  and  explosions.  Techniques  using 
spectral  ratios  have  been  applied  to  data  in  the  Western  U.S.  and  in 
Europe  with  reasonable  success,  but  the  transportability  of  the  methods 
has  not  been  satisfactorily  established.  Wiister  used  ARMA  (Auto- 
Regressive,  Moving  Average)  spectral  estimates  in  the  development  of 
discriminants.  The  discriminant  we  studied  is  based  on  autoregressive 
modeling  of  the  Lg  phase.  The  initial  data  is  from  Grant/Coyne/Ryall  CSS 
Ground  Truth  Data  Base  recorded  at  GERESS  (16  Hz  bandwidth).  Figure 

1  is  a  location  map  of  the  GERESS- Vogtiand  area.  The  data  consists  of 
explosions  and  earthquakes  in  the  Vogtiand  area  (see  Table  1  of  Appendix 

2  on  page  23  of  this  report). 

The  magnitude  range  is  1.93  -  2.15  for  the  quarry  blasts  and  1.4  -  2.37  for 
earthquakes.  The  origin  times  were  calculated  assuming  Jeffreys-Bullen 
travel  time  curves.  From  the  waveforms,  those  recorded  by  GEC2  (the  most 
reliable  station  in  the  GERESS  array)  were  used. 
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AR  (ALL  POLE^  PROCESS 


x{n)=-  ^  a{k)x{n-k)+u{n) 

k=\ 


p  is  the  order  of  the  process. 


a(k)  are  the  coefficients, 
generally  complex. 


u(n)  is  the  driving  process, 
here  assumed  to  be  white 
noise. 
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Figure  1,  The  Vogtiand  area  and  the  location  of  the  GERESS 
array(after  Wuster^  1993) 


Figure  8.  --  Location  Map  of  GERESS-Vogtland  Area 
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Seismogram  -  EVENT1 


PSD  AR(3) 


Figure  9.  --  Seismogram  and  Spectra  for  Quarry  Blast  (EVENT  1) 
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Figure  10.  --  Seismogram  and  Spectra  for  Quarry  Blast  (EVENT  3) 

52 


Seismogram  -  EVENTS 
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Figure  11.  -  Seismogram  and  Spectra  for  Earthquake  (EVENT  6) 
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Seismogram  -  EVENTS 
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Figure  12.  --  Seismogram  and  Spectra  for  Earthquake  (EVENT  9) 
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Seismoaram  -  GUMZ40 


Abs.  pole 


Absolute  pole  position  as  function  of  frequency 


Frequency[Hz] 


Figure  14.  -  The  AR(3)  model  has  one  real  pole  and  two  paired  complex 
poles.  This  figure  shows  the  three  reciprocal  poles  in  the  complex  plane. 
Angular  frequency  is  from  zero  to  ±  k,  the  latter  being  the  folding 
frequency.  The  real  pole  can  only  be  at  one  of  these  two  frequencies. 

Only  the  complex  pole  for  positive  frequency  need  be  considered  in  the 
interpretation  of  the  model.  For  AR  process  to  be  stable,  the  poles  must  all 
lie  outside  the  unit  circle,  therefore  the  reciprocal  poles  must  lie  between 
the  origin  and  the  unit  circle.  In  practice,  their  position  cannot  be  zero  or 
one.  The  nearer  the  reciprocal  pole  to  the  unit  circle,  the  “stronger”  it 
appears  in  the  AR  spectrum. 
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Figure  15.  --  Reciprocal  pole  position  vs  frequency  for  the  analyzed 
events.  Solid  circles  are  explosions  and  stars  are  earthquakes.  The  two 
open  circles  represent  surface  and  buried  explosions  at  the  Tyrnyaus  mine 
in  the  Caucasus  Mountains  of  Southern  Russia  recorded  at  station  GUM. 
This  figure  shows  a  clear  separation  of  explosions  and  earthquakes  with  the 
latter  having  broad  AR(3)  spectra  with  “weak”  poles  above  6  Hz  while  the 
explosions  all  show  much  “stronger”  poles  at  frequencies  less  than  5  Hz. 
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SUMMARY 
AR(3)  -  Lg  METHOD 


^  The  method  is  simple  and  robust. 


^  Because  spectral  amplitudes  or  spectral  ratios  are  not  used,  extensive 
calibrations  are  not  required. 


^  Only  the  vertical  record  of  Lg  is  needed.  The  method  does  not 
require  good  signal-to-noise  ratio  for  the  body  wave  arrivals. 


^  The  method  has  provided  complete  discrimination  for  all  events  we 
have  analyzed  which  were  recorded  with  16  Hz  bandwidth  and  for 
which  the  source  type  is  known. 
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APPENDIX  5.  -  INSTALLATION  OF  POSTHOLE  54000  SEISMOMETER 


Eugene  Herrin 

1.  Hole  was  prepared  using  the  technique  described  for  GS-13  installations. 
Casing  IDE  nominal  8  in.  Depth  nominal  20  ft. 

2.  54,000  PH  package  was  equipped  with  flat  nose  plate  and  standard 
armored  cable. 

3.  Casing  was  marked  N-S.  54,000  PH  was  marked  N-S. 

4.  Very  fine,  clean  sand  (sand  blasting  material)  was  introduced  into  the 
hole  to  a  thickness  of  2  to  4  inches.  Sand  was  introduced  through  a 
plastic  pipe  24  feet  long  with  a  funnel  on  top  in  order  to  minimize  dust 
in  hole. 

5.  With  the  masses  locked  the  54,000  PH  package  was  lowered  into  hole, 
positioned  approximately  in  center  of  hole  and  rotated  so  N-S  sensor 
was  lined  up  with  N-S  direction  marked  on  casing. 

6.  Sand  was  carefully  introduced  into  hole  using  plastic  pipe  until  sand 
level  was  about  1/2  way  up  on  the  54,000  PH  package. 

7.  A  steel  pipe  with  an  eccentric  on  the  lower  end  and  a  horizontal  rod  on 
the  upper  end  was  used  to  tilt  the  upper  end  of  the  54,000  PH  package 
until  the  level  bubble  on  top  of  the  package  was  centered.  Bubble  was 
observed  using  a  Nikon  spotting  scope  on  a  tripod  and  a  500,000  CP 
battery  powered  hand-lamp. 

8.  The  54,000  PH  control  box  was  connected  to  the  up-hole  end  of  the 
armored  cable.  Power  was  provided  to  the  control  box  (110  vac)  and  the 
three  masses  were  unlocked. 

9.  The  mass  position  (DC)  output  voltage  was  read  and  it  was  determined 
that  the  mass  was  captured  by  the  feed-back  process. 

10.  The  signal  output  of  the  vertical  component  was  observed  to  determine 
that  the  Z-component  was  operational. 

11.  Two  meters  were  used  to  read  the  mass  position  (MP)  output  voltage  of 
the  N-S  and  E-W  components.  In  general,  both  voltages  would  be  full- 
scale  (power  supply  voltage  about  +  13  vdc)  indicating  that  the  masses 
were  not  captured.  Polarity  was  N  and  E  positive. 


12.  The  pipe- tool  was  used  to  move  the  top  of  the  instrument  along  the  N- 
S  axis  in  the  required  direction  (toward  zero  MP  voltage).  A  steady 
pressure  was  required  to  be  applied  to  the  side  of  the  upper  end  of  the 
package  and  held  until  the  sand  could  flow  to  accommodate  a  tilt  in  the 
package.  The  process  was  repeated  until  the  N-S  mass  was  captured  (MP 
output  less  than  +  12  vdc). 

13.  This  process  can  be  time  consuming  and  frustrating.  A  great  deal  of 
time  and  pressure  was  required  to  move  the  top  of  the  package  only  a 
very  small  fraction  of  an  inch. 

14.  When  the  N-S  was  captured  the  E-W  MP  voltage  was  read.  Let  us 
assume  an  E-W  value  of  +13  volts  so  that  the  top  would  need  to  be 
moved  to  the  west.  Also  assume  a  MP  output  of  +  11  volts  indicating  a 
tilt  to  the  north. 

15.  The  pipe  tool  was  repositioned  to  push  the  top  of  the  package  to  the 
southwest.  Pressure  was  applied  until  the  N-S  MP  output  was  about  -11 
vdc. 

16.  The  tool  was  repositioned  to  push  the  top  of  the  package  to  the 
northwest  until  the  N-S  MP  was  about  +  11  vdc. 

17.  This  process  was  repeated  until  the  E-W  mass  was  captured.  A  near 
capture  condition  is  indicated  by  the  MP  output  changing  rapidly  as 
pressure  is  applied  to  the  top  of  the  hole  and  then  slowly  drifting  back  to 
the  stop  after  pressure  is  released. 

18.  Once  all  masses  were  captured  the  N-S  and  E-W  tilts  could  be  carefully 
trimmed  until  MP  output  on  both  horizontal  channels  were  +  a  couple 
of  volts  or  less. 

19.  The  signal  outputs  were  observed  to  insure  that  all  the  components 
were  measuring  seismic  activity. 

20.  The  hole  and  cable  were  put  into  the  final  operational  configuration.  If 
desired,  more  sand  can  be  introduced  around  the  54,000  PH  package  for 
stability,  but  after  the  experience  of  trying  to  tilt  the  package  using  the 
pipe  tool,  you  will  probably  conclude  that  no  more  sand  is  needed. 
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